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We present recent results obtained by the BRAHMS experiment at the Relativistic 
Heavy Ion Collider (RHIC) for the systems of Au + Au and Cu + Cu at ^/snn = 
200 GeV and at 62.4 GeV, and p + p at = 200 GeV. BRAHMS explores reaction 

dynamics and properties of the hot and high energy density matter produced in ultra- 
relativistic heavy-ion collisions versus its longitudinal expansion. Overall charged particle 
production, particle spectra over a large rapidity interval and px range are presented. 
Nuclear modification factors for Au + Au and Cu + Cu collisions are discussed. The 
observed number of charged particles produced per unit of rapidity in the central rapidity 
region indicates that a high energy density systems are produced at the initial stage of the 
Au + Au reaction. Analysis of anti-particle to particle ratios as a function of rapidity and 
collision energy reveal that particle populations at the chemical freeze-out stage for heavy- 
ion reactions at and above SPS energies are controlled by the baryon chemical potential. 
From the particle spectra we deduce significant radial expansion (/9 ^ 0.75), as expected 
for systems created with a large initial energy density. We also measure the elliptic flow 
parameter V2 versus rapidity and pt- A weak dependence on rapidity of the px differential 
V2 is observed. We present rapidity dependent p/n ratios within < y < 3 for Au + Au 
and Cu -|- Cu at ^/sNN = 200 GeV. The ratios are enhanced in nucleus-nucleus collisions 
as compared to p + p collisions. The particle ratios are discussed in terms of their system 
size and rapidity dependence. We compare Raa for Au + Au at ^/snn = 200 GeV and 
at 62.4 GeV, and for Au + Au and Cu + Cu at ^/snn = 62.4 GeV. Raa is found to 
increase with decreasing collision energy, decreasing system size, and when going towards 
more peripheral collisions. However, Raa shows only a very weak dependence on rapidity 
(for < y < 3.2), both for pions and protons. 



1. INTRODUCTION 

Reactions between heavy nuclei provide a unique opportunity to produce and study 
nuclear (hadronic) matter far from its ground state, at high densities and temperatures. 
From the onset of the formulation of the quark model and the first understanding of the 
nature of the binding and confining potential between quarks about 30 years ago, it has 
been realized that at very high density and temperature, hadronic matter may undergo 
a transition to a more primordial form of matter. This proposed state of matter named 
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the quark gluon plasma (QGP) [ 1., is characterized by a strongly reduced interaction 
among its constituents, quarks and gluons, such that the partons would exist in a nearly 
free state [ Z. Experimental attempts to create the QGP in the laboratory by colliding 
heavy nuclei have been carried out for more than 20 years. During this period, center of 
mass energies per pair of colliding nucleons have risen steadily from the ^/sNN = 1 GeV 
domain of the Bevalac at LBNL, to energies of ^/sNN = 5 GeV at the AGS at BNL, and 
to ^/SNN = 17 GeV at the SPS at CERN. Although a number of signals suggesting the 
formation of a very dense state of matter have been observed, no solid evidence for QGP 
formation was found at these lower energies. 

In mid- August 2001 systematic data collecting by the four RHIC experiments, namely 
BRAHMS [01, PHENIX [g], PHOBOS [ | and STAR [i], began at the energy of 
^/SNN = 200 GeV. The RHIC operations started a new era of studies of ultra-relativistic 
nucleus-nucleus collisions. One of the main results obtained up to the present at RHIC is 
the observation of significant elliptic flow in central Au + Au collisions. The large flow 
signal, that is consistent with the hydrodynamic evolution of a perfect fluid [I71IHI 1^ ITU], 
indicates a strongly interacting QGP, contrary to initial expectations. Also observed is 
strong jet suppression [E], which is predicted within the Quantum Chromo-Dynamics 
theory (QCD) as a consequence of the creation of a dense colored medium [[I21- This 
last observation was supplemented by d + Au measurements showing the absence of 
suppression, but rather a Cronin type enhancement at the central rapidity region, thus 
excluding an alternative interpretation of the suppression in terms of initial state parton 
saturation (CGC) effects. Another measurement that strongly supported the scenario of 
jet quenching was the disappereance of the away-side jet [ UHl Ej in central Au -|- Au 
collisions, whereas for d -|- Au and peripheral Au + Au near side and away side back-to- 
back jet correlations have been measured. Its large rapidity range and px coverage allows 
BRAHMS to study the properties of the produced medium as a function of its longitudinal 
expansion. The measurement of rapidity evolution of the nuclear modification factor for 
d + Au performed by BRAHMS shown that at more forward rapidities the hadronic yields 
are suppressed as compared to scaled p + p interactions [HH]- The suppression was even 
stronger for central collisions. Both observations can be quantitatively described within 
the framework of the CGC [US CHI- 

2. BRAHMS DETECTOR SETUP 

The BRAHMS (Broad RAnge Hadron Magnetic Spectrometers) [El? experimental 
setup consists of a set of global detectors and two spectrometer arms: a Mid-Rapidity 
Spectrometer (MRS) that operates in the polar angle range from 30 to 90 degrees, and a 
Forward Spectrometer (FS) that operates in the range between 2.3 and 30 degrees. Global 
detectors are used to measure the global features of the collision such as the overall particle 
multiplicity and collision centrality, position of the collision vertex and, more recently, in- 
formation on the reaction plane orientation which can be used for azimuthal flow analysis. 
For the momentum measurements in the MRS we use two tracking devices and one dipole 
magnet. Particle identification (PID) is done by Time-of-Flight (TOF) measurement and 
by using threshold Cherenkov detector C4. In the FS two Time Projection Chambers 
(TPC) and three Drift Chambers (DC) deliver particle track segments to allow excellent 
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momentum resolution using three dipole magnets. PID in the FS is provided by TOF 
measurements for low (TOFl) and medium (T0F2) particle momenta. High momentum 
particles are identified using a Ring Imaging Cherenkov detector (RICH). The BRAHMS 
PID ability is summarized in Tabled 

Table 1 



Upper range of the momentum for 2cr separation (in GeV/c) 



< < 1.0 


1.5 <r] < 4.0 




TOFW T0FW2 


C4 


TOFl T0F2 


RICH 


K/n 
K/p 


2.0 2.5 
3.5 4.0 


9.0 


3.0 4.5 
5.5 7.5 


25.0 
35.0 



3. OVERALL BULK CHARACTERISTICS 

The multiplicity distribution of emitted particles is a fundamental observable in ultra- 
relativistic collisions. It is sensitive to all stages of the reaction and can address issues such 
as the role of hard scatterings between partons and the interaction of these partons in the 
high-density medium [UHlIiniEn!- Figure [T] shows the measured pseudo-rapidity density of 
charged hadrons, dNch/d?7 , over a wide range of rj for central (0 — 5%) Au + Au collisions 
at = 200 GeV, 130 GeV and 62.4 GeV. For central collisions at = 200 GeV 

we observe about 4500 charged particles within the rapidity range covered by the detection 
system and dNch/drj |^=o = 625 ± 56. The latter value exceeds the particle production 
per participant pair observed in elementary p + p collisions at the same energy by 40 - 
50% [1^. This means that nucleus- nucleus collisions at the considered energies are far 
from being the simple superposition of elementary nucleon-nucleon collisions. 

The measurement of charged particle density dNch/dr] can be used to estimate the 
so-called Bjorken energy density, e The formula 



_3 {Et) dN,h .... 
^ " 2 vri?2r„ dr^ 

provides the value of about 4 GeV/fm'^. To obtain this result we assumed that Tq = 1 fm/c, 
[Et) = 0.5 GeV and R = 6 fm. The factor 3/2 is due to the assumption that the 
charged particles carry out of the reaction zone only a fraction (2/3) of the total available 
energy. The more refined results obtained from identified particle abundances and particle 
spectra leads to somewhat larger values of e, namely 5 GeV/fm^ at ^/sNN = 200 GeV, 
4.4 GeV/fm3 at = 130 GeV, and 3.7 GeV/fm^ at = 62.4 GeV, [EHj. All 

of these values significantly exceed the predicted energy density e ^ 1 GeV/fm^ for the 
boundary between hadronic and partonic phases of nuclear matter [1^. 

3.1. Hadrochemistry with BRAHMS data 

BRAHMS measures anti-particle to particle ratios for pions, Nt^-/Nt^+, kaons, N^- /Nx+, 
and protons, Np/Np, over a large rapidity interval. The results, for Au -|- Au collision at 
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Figure 1. Distributions of dNch/drj measured by BRAHMS for — 5% central Au + Au 
reactions, at ^/sNN = 200 GeV (solid circles), ^/sNN = 130 GeV (open squares) and 
y/SNN = 62.4 GeV (solid squares). Only statistical errors (usually smaller than the 
symbol size) are shown. 



^/sNN = 200 GeV and ^/snn = 62.4 GeV are presented in Figure El Whereas N^,- /Nt,+ 
stays constant and is equal to 1 over the covered rapidity range (0 < y < 3.2) for both 
energies, the Nk-/Nk+ and Np/Np ratios drop significantly with increasing rapidity. For 
y/s^jy = 200 GeV the Nk-/Nk+ and Np/Np ratios are equal, respectively to 0.95 and 
0.76 at ?/ ~ 0, and reach values of 0.6 for Nk-/Nk+ and 0.3 for Np/Np around rapidity 
3. Comparing ^/snn = 200 GeV and y/s^N = 62.4 GeV we observe a 11% and 40% 
decrease of the ratio at mid-rapidity and a 43% and 93% decrease at j/ ~ 3 for kaons and 
protons, respectively. The large change in the Np/Np ratio from ^/sNN = 200 GeV to 
^/sNN = 62.4 GeV is due to the fact that the rapidity ?/ ^ 3 for the lower energy cor- 
responding to the fragmentation region, whereas for ^/sNN = 200 GeV, y ~ 3 is located 
about 1 unit of rapidity below the maximum in the net-baryon distribution [I25j. 

Figure ISl shows the Nk-/Nk+ ratio as a function of the corresponding Np/Np for var- 
ious rapidities. The presented results were obtained for central collisions at three RHIC 
collision energies. The AGS and SPS results are plotted for comparison. There is a 
striking correlation between the RHIC/BRAHMS kaon and proton ratios over 3 units 
of rapidity. It is worth noting that the BRAHMS forward rapidity data measured at 
^/sNN = 62.4 GeV overlap with the SPS points that were measured at much lower energy 
but at mid-rapidity. The solid line in Figure 3 shows a fit with a statistical model to the 
^/sNN = 200 GeV results, only, assuming that the temperature at the chemical freeze-out 
is 170 MeV [1211 12Z]- It is seen that the data are very well described by the statistical 
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Figure 2. Ratios of anti-particles to particles (pions, kaons and protons) as a function of 
rapidity for ^/snn = 62.4 GeV, ^/sNN = 200 GeV. Statistical and systematic errors are 
indicated. 



model over a broad rapidity range with the baryon chemical potential changing from 27 
MeV at mid-rapidity to 140 MeV at the most forward rapidities. Using simple statistical 
models at the quark level with chemical and thermal equilibrium, the ratios can be written 



Np 



-6Mu,d/r 



K- 



Nk+ 



-2(M„,d-M»)/T 



(2) 



where and T are the chemical potential and temperature, respectively. Substituting 
= into eqs. 0, one gets Nk-/Nk+ = [N^/N^'^/^ . This relation, represented by 
the dotted line on Figure 121 does not reproduce the observed correlation. The data are, 
however, well fitted by the function N^- /Nx+ = [iVp/iVp]^/'^ (dashed line) which can be 
derived from eqs. (j2I) assuming /z^ = l/4/x„ rf. 

Recently, STAR and NA49 have measured mid-rapidity ratios A/A, S/S and VL/VL 
versus Np/Np for a set of energies from ^/sNN = 10 GeV up to ^/snn = 200 GeV. These 
preliminary results can also be well described within a statistical model of chemical and 
thermal equilibrium at the quark level and confirm the strong correlation between fis and 
^u,d derived from BRAHMS data. 

3.2. Radial flow 

The properties of matter in the latest stage of the collision when the interactions be- 
tween particles cease (kinetic freeze-out) can be studied from the shapes of emitted particle 
spectra. These shapes depend in general on the temperature of the emitting source and 
on the collective flow. For central collisions where one should not expect any azimuthal 
dependence, only the so-called transverse flow is important [E]- In the blast-wave ap- 
proach [I2H1 the spectrum is parametrized by a function which depends on the freeze-out 
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Figure 3. Correlation between NK~/Nj^+arLd Np/Np. The solid curve refer to statistical 
model calculation with a chemical freeze-out temperature of 170 MeV. 



temperature, Tjo, and on the transverse expansion velocity, Pt- Figure El shows results 
obtained from the analysis of the Au + Au reactions at ^/sNN = 200 GeV. The blast- 
wave model was used to fit, simultaneously, the vr"*", vr", K"*", K~, protons and anti-protons 
spectra. The results are plotted versus the number of participants, Npart- 

One can see that the kinetic freeze-out temperature decreases with centrality from about 
140 MeV for the 40 — 50% centrality bin to about 120 MeV for the most central collisions. 
The latter value is lower than the temperature of chemical freeze-out indicating that, as 
expected, the freeze-out of particle ratios occurs earlier than the kinetic freeze-out. The 
reversed trend is seen for the expansion velocity which is equal to about 0.65c and 0.75c 
for 40 — 50% and — 5% collision centrality, respectively. For comparison we show also 
the result for the central Au + Au collisions at a/s^vtv = 62.4 GeV. Although for the same 
value of Npart the Tfo value for both energies is very similar, a reduction by about 20% 
in the expansion velocity is observed. 

Figure El shows the dependence of Tfo and Pt as a function of rapidity. The observed 
variations of Tfo and Pt on Npart and rapidity are consistent with the hydrodynamic 
description in which the radial flow is the result of outwards gradients of pressure that 
exist in the expanding matter during the whole evolution. Thus the speed of expansion 
should increase with the density of the initially created system. The transverse expansion 
velocity is larger than that observed at SPS energies, which is consistent with a large 
initial density of the system created at RHIC. 
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Figure 4. Kinetic freeze-out temperature 
and transverse flow velocity at mid-rapidity 
as a function of centrality for Au + Au at 
y/SNN = 200 GeV (dots). For comparison we 
show the result for — 10% central Au + Au 
collisions at ^/snn = 62.4 GeV (stars). 



Figure 5. Kinetic freeze-out temper- 
ature and transverse flow velocity for 
central Au + Au collisions at ^/sNN = 
200 GeV function of rapidity. 



3.3. Elliptic Flow 

A powerful tool for studying the dynamics that drive the initial evolution of systems 
created in heavy ion reactions is the analysis of the azimuthal distribution of the emitted 
particles relative to the reaction plane [E]. The triple differential distribution of emitted 
particles can be factorized as follow 

^ [I + 2vi{y,pT)cos{(f)-(f)r)+2v2iy,PT)cos2{(j)-(j)r) + ....), (3) 



dydprdcj) dydpx 2n 

where and (pr denote the azimuthal angles of the particle and of the reaction plane, 
respectively. The first factor depends only on y and px and the second factor represents 
the expansion of the azimuthal dependence into a Fourier series. The coefficients of the 
first (fi) and the second (f2) harmonics are called direct and elliptic fiow parameters, 
respectively, and in general they are functions of rapidity and transverse momentum. 
Calculations based on hydrodynamic models [HHI show that elliptic fiow is substantially 
generated only during the highest density phase, before the initial spatial anisotropy 
of the created medium disappears. Thus V2 is sensitive to the very early phase in the 
system evolution and relatively insensitive to the late stage characterized by the dissipative 
expansion of the hadronic gas. 

Elliptic fiow has been extensively measured by the STAR, PHENIX and PHOBOS 
experiments [HUHlini. Generally, V2 is an increasing function of px up to 1.5 GeV/c, at 
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Figure 6. Pion elliptic flow parameter pt dependence measured for Au + Au at ^/sNN 
200 GeV at r/ = and at r/ = 3.4. 



which point it saturates. Up to roughly 1.5 GeV/c in p^, hydrodynamic calculations show 
good agreement with experimental data for the f 2 dependence on px and centrality. 

BRAHMS has measured the pt dependence of V2 at a number of rapidities [1211, 
results for pions at 77 = and at = 3.4 are shown in Fig. IHl The elliptic flow parameter, 
V2, for identified pions is an increasing function of pt at both rapidities. The dependence 
on pseudo-rapidity is very small. These results are very similar to those obtained for 
charged hadrons [12^. 



4. BARYON TO MESON RATIOS 



With its excellent particle identification capabilities, BRAHMS can study the pt and y 
dependence of hadron production. Preliminary results [EOlEIl indicate that for Au + Au 
reactions in the intermediate pt region the proton to meson ratio is significantly higher 
that one would expect from the parton fragmentation in vacuum process. Several theoret- 
ical models have been proposed to explain the observed enhancement. These range from 
models that explore the partonic interactions of quark hadronization and quark coales- 
cence [021 ESI Ell to models that incorporate novel baryon dynamics [EniEHl ■ The recent 
experimental data obtained for the p + p, Au -|- Au, and Cu + Cu colliding systems are 
expected to result in a better understanding of the underlying physics. Fig. [7| shows a 
recent investigation by BRAHMS [E3 of P to vr" ratios at mid-rapidity (circles) and at 
pseudo-rapidity r] ^ 3.2 (squares) for Au + Au collisions at ^/snn = 200 GeV for the dif- 
ferent centrality classes indicated on the plot. The data reveal a smooth increase of p/vr" 
ratios from peripheral to central collisions, however, the centrality dependence is stronger 
at mid-rapidity than at forward rapidity. The maxima in the p/vr" ratios are lower at 
forward rapidity as compared to mid-rapidity. Figure El shows the p/tt" centrality depen- 
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Figure 7. The p to vr^ ratios at 
mid-rapidity (circles) and at pseudo- 
rapidity r] ^ 3.2 (squares) for differ- 
ent centralities of Au + Au collisions 
at ^/s^ = 200 GeV. 



Figure 8. The averaged p/vr~ versus (Npart) for 
Au + Au (open symbols) and Cu + Cu (solid 
symbols) at ^/sNN = 200 GeV, for ?/ ^ (cir- 
cles) and and for r] ^ 3.2 (squares). 



dence for Au + Au (open symbols) and Cu + Cu (solid symbols) at y/s^N = 200 GeV. 
The data ior y = and rj ^ 3.2 are represented by circles and squares, respectively. 
One can see the strong increase of the p/vr^ ratios as a function of Np^^ in the range 
< Npart < 60. For Npart > 60 the dependence starts to saturate and the ratios reach 
values of about 0.6 and about 0.25 for central collisions at mid- and forward rapidities, 
respectively. For peripheral Au + Au colhsions the data approach the p -l- p results. It is 
important to note that the Au -|- Au and Cu + Cu ratios are consistent with each other 
when plotted versus (Npart), indicating that the p/ir^ ratio is controlled by the initial size 
of the created systems. Fig. |H1 shows the comparison of BRAHMS and PHENIX data 
for the ratio of protons to positive pions measured at y = 0. The parton coalescence [ 
and recombination models describe the observed ratios well at mid-rapidity, as 
shown in the figure by the dashed curve. The three-dimensional hydrodynamic model 
[ IHS] (solid curve) cannot reproduce the observed shape of the p/vr^ ratio, however it 
reproduces the overall level of enhancement. The current hydrodynamic calculations do 
not show a rapidity dependence which is contrary to the experimental observation (see 
Fig. El). This is presumably because the hydrodynamic models assume constant baryon 
chemical potential versus rapidity. 
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Figure 9. The p/ir'^ ratio versus px for Au + Au collisions at ^/sN^ = 200 GeV measured 
at mid-rapidity by BRAHMS (circles) and by PHENIX (triangles). Squares represent 
BRAHMS results at ~ 3.2. The plotted lines show model calculations (see text). 



5. HIGH pt suppression 

Particles with high pt (above 2 GeV/c) are primarily produced in hard scattering pro- 
cesses early in the collision. In high energy nucleon-nucleon reactions hard scattered 
partons fragment into jets of hadrons. However, in nucleus-nucleus collisions hard scat- 
tered partons might travel in the medium. It was predicted that if the medium is a QGP, 
the partons will lose a large fraction of their energy by induced gluon radiation, effec- 
tively suppressing the jet production [EH]. Experimentally this phenomenon, known as 
jet quenching, will be observed as a depletion of the high px region in hadron spectra. 

The measure most commonly used to study the medium effects is called the nuclear 
modification factor, Raa- It is defined as the ratio of the particle yield produced in 
nucleus-nucleus collision, scaled with the number of binary collisions (Ncoii), and the 
particle yield produced in elementary nucleon-nucleon collisions: 

^ YieldjAA) 
N,,uxYteld{NN)- 

At high Pt the particles are predominantly produced by hard scatterings and in the 
absence of nuclear effects (when the nucleus-nucleus collision reduces to the superposition 
of elementary collisions) we expect Raa to be 1. At low pr, where the production rate 
scales rather with Np^rt, Raa should converge to '^part/'^coii which is roughly 1/3 for 
central Au + Au collisions at the top RHIC energy. Raa < 1 at high pt will indicate 
the suppression which, as has been discussed, can be attributed to the jet quenching 
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Figure 10. Raa measured at r/ = and rj = 1 ioi Au + Au at ^/snn = 200 GeV (upper 
row) and at ^/snn = 62.4 GeV (bottom row), for different centrality bins indicated on 
the plot (p + p reference for -v/sjvtv = 62.4 GeV is based on ISR collider data [ 40j). 



phenomenon. At SPS there is no suppression, in fact it is well known that there is 
enhancement for > 2 GeV/c. This so-called Cronin effect is attributed to initial 
multiple scattering of reacting partons. 

Another variable used to quantify nuclear effects, one that does not depend on the 
elementary reference spectra, is Rcp, defined as the ratio of Raa for central nucleus- 
nucleus collisions to Raa for peripheral collisions. The idea of using Rcp is based on the 
expectation that any nuclear modification in peripheral collisions will be insignificant. We 
will show that the latter statement is not true for Rcp = -RauA«(0— 10%)/i?^„^„(40— 60%). 

5.1. Raa evolution on collision centrality and collision energy for Au -|- Au 
and Cu -|- Cu systems 

The large set of data collected during the RHIC run 4 has allowed us to carry out the 
study of the Raa evolution on collision centrality and collision energy for two colliding 
systems, namely Au + Au and Cu + Cu [HI]. Figure ITUl shows Raa measured at = and 
?7 = 1 for charged hadrons produced in Au -|- Au reactions at ^/snn = 200 GeV (upper 
row) and at ^/sJfN = 62.4 GeV (bottom row), for different collision centralities indicated 
on the plot. For the most central reactions Raa shows suppression for both energies, 
however, the suppression is significantly stronger at the higher energy. We observe a 
smooth increase of Raa towards less central collisions, for ^/snn = 200 GeV, resulting in 
approximate scaling with Neon for pt > '2 GeV/c for the 40 — 50% centrality bin. However, 
at y/SNN = 62.4 GeV, Raa ~ 1 for more central collisions, and a Cronin peak is clearly 
visible already for the 20 — 40% centrality class, where Raa reaches value of about 1.3 
in the pt range between 2.0 and 3.0 GeV/c. These observations at ^/snn = 62.4 GeV 
are qualitatively consistent with a picture in which there are two competing mechanisms 
that influence the nuclear modification in the intermediate and high pt range, namely: 
jet quenching that dominates at central collisions and Cronin type enhancement (/cy 
broadening or/and quark recombination) that prevails for the more peripheral collisions. 
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Figure 11. Raa measured at mid-rapidity for Au + Au (upper row) and for Cu + Cu 
(bottom row) at ^s^n = 62.4 GeV for different centrality ranges indicated on the plot. 



Figure ITD presents similar comparison as shown on Fig. ^Jbut this time we compare 
two different systems, namely Au + Au and Cu -|- Cu colliding at the same energy 
(a/s/vtv = 62.4 GeV). For Cu + Cu at y/s^N = 62.4 GeV the same trend of increasing Raa 
with decreasing collision centrality is seen. However, now the Cronin type enhancement 
is present already for the most central collisions. 

Summarizing the whole set of observation we conclude that the level of suppression of 
the inclusive hadron spectra produced in nucleus-nucleus collisions at RHIC energies in 
the pt range above 2 GeV/c increases with increasing collision energy, collision centrality 
and with the size of the colliding nuclei. The dependency on the last two variables can 
be replaced by only one dependency on 'Npart [1^. 

5.2. Raa for identified hadrons at forward rapidity for Au + Au at y^J/vw = 
200 GeV 

The BRAHMS collaboration has already noted that the nuclear modification factor in 
central Au + Au at ^/snn = 200 GeV at r/ = 2.2 is comparable to that measured at 
mid-rapidity [ |Ilj for the same system and energy. However, the data did not clearly 
identify the mechanism responsible for the observed effect. In this section we present 
more exclusive analysis of nuclear effects at forward rapidity for identified pions, kaons 
and protons. 

Figure [T21 shows Raa for identified hadrons at rapidity y ~ 3.2 for — 10% central 
Au + Au events at ^/snn = 200 GeV. The shaded band around unity indicates system- 
atic error associated with the uncertainty in the number of binary collisions. Both Raa 
and Rcp show suppression for pions (left panel) and for kaons (middle panel), however 
for protons (right panel) Rcp shows suppression whereas Raa reveals a Cronin type en- 
hancement with the peak at ~ 2 GeV/c. The difference between Raa and Rcp for 
protons is striking, indicating significant enhancement of protons (in respect to p + p) 
for the 40 — 60% collision centrality used in the definition of Rcp- The same misleading 
behavior of Rcp is seen for charged hadrons when comparing evolution of Raa and Rcp 
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Figure 12. BRAHMS Raa (upper row) and Rcp (bottom row) for pions, kaons and 
protons, measured y ~ 3.2 in Au + Au ^/snn = 200 GeV. The p + p reference was also 
measured by BRAHMS (see 



on T] for Au + Au at y/sj^ = 200 GeV and y/sj^ = 62.4 GeV [liT]. 

Figure IT^ shows the nuclear modification factors calculated for (vr"*" + 7r~)/2 (left panel) 
and (p + p)/2 (right panel), respectively, at ?/ ~ 3.2, for central Au + Au reaction. For 
the comparison we also show the Raa values measured by the PHENIX Collaboration 
at mid-rapidity. The Raa measured for pions shows strong suppression (by factor of 
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Figure 13. Comparison of Raa measured for central Au + Au collisions at ^/snn 
200 GeV, at mid-rapidity and y ^ 3 for pions (left panel) and protons (right panel). 



about 3 for 2 < < 3 GeV/c), both at mid- and at forward rapidity. The consistency 
between mid- and forward rapidity is seen also for protons, but in this case, Raa reveals 
a Cronin peak around pt = 2 GeV/c. The similarity between Raa at mid- and forward 
rapidity observed simultaneously for pions and protons suggests the same mechanisms is 
responsible for the nuclear modifications within the studied rapidity interval. 

It has been predicted [HHl that the magnitude of jet quenching should depend on both 
the size and the density of the created absorbing medium, thus it is interesting to study 
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the dependence of Raa on centrality. In Figure [T31 we plot the averaged Raa measured 
for pions versus {Npart) for mid-rapidity (sohd circles) and for forward rapidity (open 
squares). The averaging was performed in the pt range from 2 GeV/c to 3 GeV/c. Once 
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Figure 14. Averaged Raa m the range 2.0 < px < 3.0 GeV/c at mid-rapidity (PHENIX) 
and at forward rapidity as a function of collision centrality, expressed by the number of 
participants. 



again, the mid- and forward rapidity pion suppression for the most central Au + Au 
reaction are found of the same strength. However, the Raa measured at forward rapidity 
shows significantly stronger rise towards peripheral collisions as compared to Raa at mid- 
rapidity, differing on the level of 35% for (Npart) ~ 100. The trend is consistent with the 
model of parton energy loss in a strongly absorbing medium [031 US]- In this picture, for 
y = 0, the jet emission is dominated by the emission from the surface which quenches the 
dependence of Raa on the system size. On the other hand, for y ~ 3, the transition from 
surface to volume emission can occur, which leads to a stronger dependence on Np^rt- 

6. SUMMARY 

The results from BRAHMS and the other RHIC experiments clearly show that studies of 
high energy nucleus-nucleus collisions have moved to a qualitatively new physics domain. 
The collisions are characterized by a high degree of reaction transparency leading to the 
formation of a near-baryon-free central region with approximate balance between matter 
and antimatter. From the measurement of charged particle multiplicities in this region 
lower limits for the energy density at Tq = 1 fm/c have been determined as 5 GeV/fm^ and 
3.7 GeV/fm^, for central Au -|- Au reactions at ^/sNN = 200 GeV and ^/sNN = 62.4 GeV, 
respectively. Therefore the conditions necessary for the formation of a deconfined system 
appear to be well fulfilled at RHIC energies. Analysis within the statistical model of the 
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relative abundances of K~, K^, p and p suggests equilibrium at a chemical freeze-out 
temperature of 170 MeV, with a noticeably strong correlation between the strange quark 
and baryon chemical potentials. Analysis of particle spectra within the blast-wave model 
indicates a kinetic freeze-out temperature on the order of 120 MeV and a large transverse 
expansion velocity, which is consistent with the high initial energy density. 

The measurement of the elliptic flow parameter f 2 versus rapidity and pt shows weak 
dependence of the pt differential V2 on rapidity. The p/ir ratios were measured within 
< 77 < 3 for Au + Au and Cu + Cu at = 200 GeV and ^saw = 62.4 GeV. The 

obtained results reveal strong enhancement of baryon to meson ratios in nucleus-nucleus 
collisions as compare to p + p collisions. Models that incorporate an interplay between 
soft and hard processes can describe the data at mid-rapidity. We compared Raa for 
Au -I- Au at y/SNN = 200 GeV and y/SNN = 62.4 GeV, and for Au -|- Au and Cu -|- Cu 
at y/SNN = 200 GeV. The general observed trend is that Raa increases with: decreasing 
collision energy, decreasing system size, and when going towards the more peripheral 
collisions. For Au -|- Au central collisions at ^/snn = 200 GeV, Raa shows very weak 
dependence on rapidity (in < y < 3.2 interval), both for pions and protons. 
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